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Chapter 1 
Introduction 


An analysis is developed for predicting the electromagnetic fields (EM) 
which are scattered from a 3-dimensional (3-D) perfectly conducting S- 
shaped open-ended semi-infinite waveguide cavity with a planar interior 
termination when it is illuminated by an external plane wave. The analysis 
is based on a hybrid asymptotic modal approach which has been discussed 
extensively in [1,2, 3, 4], 

Briefly, this hybrid asymptotic-modal, or more simply the hybrid modal 
approach, models the S-shaped rectangular cavity by smoothly joining to- 
gether straight and annular rectangular waveguide sections as shown in Fig- 
ure 1.1. The scattering matrices which describe the reflection and transmis- 
sion coefficients of the waveguide modes reflected and transmitted at each 
junction between the different waveguide sections, as well at the scattering 
from the edges at the open end are found via asymptotic high frequency 
approximations such as the geometrical and physical theories of diffrac- 
tion used in conjunction with the equivalent current method. The multiple 
wave interactions between the various waveguide sections can be accounted 
for in a straightforward fashion, in terms of these junction scattering (re- 
flection and transmission) matrices, via the self-consistent multiple scat- 
tering method [1,2, 3, 4]. It is noted that the asymptotic high frequency 
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(a) oblique view 


JUNCTION (J) 



Figure 1.1: S-shaped cavity with rectangular cross-section formed by joining 
straight and annular waveguide sections. 






based approach is far more efficient for calculating the modal reflection 
and transmission coefficient matrices for the waveguide junctions than is 
the classical mode matching technique, as the latter requires a numerical 
matrix inversion making it cumbersome and inefficient. Furthermore, these 
modal coefficients obtained via asymptotic high frequency approximations 
generally have a relatively simple form thereby lending physical insight into 
the modal reflection and transmission process. It is the hybrid combination 
of high frequency and modal techniques that makes the present approach 
a hybrid asymptotic-modal approach. 

The modes in the straight and annular rectangular waveguide sections 
are known analytically, and the decomposition of these modal fields into 
their corresponding equivalent sets of modal rays can be performed easily. 
However, consistent with some previous experience [1], it is seen that the re- 
flection coefficient matrix elements are negligible in comparison to the trans- 
mission coefficient matrix elements characterizing the smooth junctions be- 
tween the straight and annular waveguide sections. Hence, it is found that 
the modal reflection effects can be ignored in this analysis thereby requiring 
one to essentially retain only the modal transmission matrix. 

It is noted that in this multiple scattering method, the total scattered 
field consists of a superposition of three contributions, namely: 

(i) contribution from the diffraction of the incident wave by the edge at 
the open end; 

(ii) contribution from the scattering by the interior termination, which 
arises from part of the incident energy which couples via the open end into 
the S-duct (or cavity) and then reflects from the termination to partly ra- 
diate out through the open end. The remaining energy is scattered back 
into the cavity to again reflect from the termination, and so on. Since 
the reflection from the planar perfectly-conducting termination is generally 
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more significant than the reflection back into the cavity from the open end, 
the multiple wave interactions between the termination and the open end 
may be ignored. Also, if the reflection from the termination can be made 
small by making it an appropriate planar dielectric termination (say), then 
the interactions between the open end and the low reflection type termina- 
tion can again be ignored! While the reflections of waves at each junction 
between the waveguide sections can contribute to the net scattering from 
the interior, this effect is negligible in comparison to the other effects as 
mentioned earlier, and it is therefore ignored making the only contribution 
to the scattering from the interior of the cavity as being due to the termi- 
nation. Nevertheless, it is noted that any of these effects which have been 
ignored here can be put back directly via the multiple scattering method 

[1.2. 3. 4] , if it is so desired; 

(iii) contribution from the scattering by any other external features of 
the S-duct. This effect is not of concern to the present work and it will 
therefore not be included. 

Since the diffraction from the open end and the modal reflection from 
the interior planar termination have already been found for the rectangu- 
lar waveguide case, and the corresponding scattering matrices have been 
given explicitly for these two junctions in [2,4], it only remains to find the 
modal transmission coefficient scattering matrix for the junction between 
a straight and annular rectangular waveguide section. Before finding the 
latter, the general formulation for the total scattered field in terms of the 
multiple scattering method is first given in Chapter 2. This multiple scat- 
tering method directly employs the scattering matrix for the open end, the 
transmission matrices for the junction between straight and annular waveg- 
uide sections, and the reflection matrix for the planar interior termination 

[1.2. 3. 4] . The details involved in the development of the modal transmis- 
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sion (scattering) matrix for the junction between a straight and annular 
waveguide section are then presented in Chapter 3. Finally, the results for 
the open end and the termination scattering matrices are combined with 
the junction transmission scattering matrix developed in Chapter 3, in ac- 
cordance with the formulation given in Chapter 2 for the total scattered 
field; the numerical results based on this total scattered field are presented 
in Chapter 4. An e iwt time convention is assumed and suppressed in the 
following analysis. 
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Chapter 2 

Formulation Using the 
Self-Consistent Multiple 
Scattering Method 


As mentioned in Chapter 1, an incident plane wave excites the fields inside 
the duct configuration as shown in Figure 1.1. These fields can be expanded 
into the natural waveguide modes of each uniform section. Transmission 
and reflection of waves at a junction as well as propagation through a 
given waveguide section are all described in terms of reflection, transmission 
and propagation matrices. The elements of these matrices are found us- 
ing asymptotic high frequency approximations such as the geometrical and 
physical theories of diffraction in conjunction with the equivalent current 
method. The multiple wave interactions between the different junctions 
can be accounted for using the self-consistent Multiple Scattering Matrix 
(MSM) formulation [1,2, 3, 4] which makes use of these reflection, transmis- 
sion and propagation matrices. 

As a quick review of the MSM formulation, consider a a plane wave 
incident on an open ended parallel plate duct as shown in Figure 2(a). 
This 2-D scattering example is considered due to its simplicity in geometry 
and in formulation. The MSM result obtained can then be extended to 
a rectangular S-shaped duct as will be shown later in this chapter. The 
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scattering matrices involved in a rectangular duct are more complex than 
those in a parallel plate duct. Reference [2] gives a detailed discussion and 
derivations of these scattering matrices. For the 2-D scattering example, 
let the symbol U be the z-directed E-field of the TE case and the E-field 
of the TM case, respectively; this incident field can be written as 

U inc = zUince’ jk( * co '* +vtine) (2.1) 

Also, let U* refer to the field scattered into the exterior region by the waveg- 
uide cavity. Thus, U * represents the complex amplitude of the scattered 
electric field E * — zE * for the TE case and the scattered magnetic field 
H‘ = zH‘ for the TM case. Let 

v = u‘+v:+u:„ ( 2 . 2 ) 

where U * is the component of the field scattered into the exterior by just 
the edges at the open end, and U* is the contribution to the field scattered 
into the exterior through reflection from the interior termination. Also, U* xt 
is the contribution to the scattered field due to any other external features 
of the duct configuration. This contribution is not of interest here and thus 
is neglected. In the far zone of the open end of the cavity, one can write 

V, = (2.3) 

v r 

e -jkr 

u; = CWe,fc(M)^ (2.4) 

V r 

where (r, 0) are shown in Figure 2.1(a), and f e ,h and g e> h are unknown 
pattern functions yet to be found. 

From the self-consistent multiple scattering method [1,2, 3, 4], one can 
express the total far zone U‘ (minus U* xt ) in terms of the generalized scat- 
tering matrices introduced as 

[t/'(iO] = {[s„] + [s„)[P)[r B )[P]((/] - [s M )(i')[r*][P])- 1 [5«]} 

(2.5) 
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(a) parallel plate duct 
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(b) 2-dimensional S-shaped duct 
Figure 2.1: 2-dimensional parallel and s-shaped duct models 




In (2.5), the matrix [P] denotes the phase delay of the modes as they 
propagate between junctions A and B; this [P] is a diagonal matrix. The 
[Tg] in (2.5) denotes the total reflection coefficient scattering matrix at 
junction B which describes the modes reflected back into region 2 due to 
modes incident at B from region 2. 

Note that if there was no termination, then [Tb] would be a null matrix 
and the scattered field U * in (2.5) would come only from the diffraction of 
the incident wave at the front edges as described by [Sn]. Therefore, [Sn] is 
the only scattering matrix present when there are no internal discontinuities 
or terminations within the cavity. From (2. 3), (2.4) and (2.5), one obtains 
that 


U = [ 5 „] 

= (2-6) 

where the matrix [Sn] has only a single element Sn; also, it follows from 
(2.2) that 

9,,h = (s. 2 ]|PJ[r B ][P]([/)-[s I 2][i’]|r B ][P])- , !s 2 .] (2.7) 

If the duct configuration is changed from a parallel plate section as 
shown in Figure 2.1(a) to an S-shaped duct section as shown in Figure 
2.1(b), the MSM result of (2.5) is modified by simply replacing the termi- 
nation reflection scattering matrix, [P b], with an effective scattering matrix, 
[Sr]. Since the front section of an S-shaped duct is a parallel plate section, 
it would be useful to define an effective scattering matrix, [Sr], at the end of 
the first section to account for all of the rest of the duct sections including 
the final termination. The MSM results of (2.5) can then be applied to find 
the scattered field of an S-shaped duct in terms of this [Sr]. For a smoothly 
joined junction, it has been shown that the modal reflection coefficients are 
negligible compared to the modal transmission coefficients [1]. Thus the 
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scattering matrix, [Sr] at junction 1 due to sections 2, 3, 4, 5 can be written 


as: 


[Sr] = WirsPl (2.8) 

with 

(4 = W][r«][P4][r3 4 ][p 3 ][r M l[/> ! ][J’n], (2.9) 

and [A] 1 is the transpose of [A]. [A] is a matrix which traces the modal 
field structure from junction 1 to the termination. [I^b] is the reflection 
coefficient scattering matrix due to the termination, and [Pj] is the propa- 
gation phase matrix due to section i. [T^] is the transmission matrix at the 
junction between sections i and j. The subscript of [Tjj] means propagation 
is from duct section i to section j . It is also noted that [A] in (2.9) gives a 
physical insight of how modes propagate in the duct structure. (2.9) should 
be read from right to left. First, the modal fields pass through a junction 
from section 1 to 2. The change in the modal field structure is accounted 
for by the modal transmission matrix [Xi 2 ]. Next, the modes propagate 
through section 2, and their phase delay results in a [P 2 ]. A subsequent 
transmission through a junction gives another transmission matrix [T23] and 
so on, until the modal fields hit the termination. 

The MSM formulation in (2.5) and the effective scattering matrix in 
(2.8) can also be applied to a 3-dimensional S-shaped duct. In Reference 
[2], an analysis has been developed for predicting the near and far zone fields 
scattered by a rectangular waveguide terminated with a planar impedance 
surface. However, if one uses the effective scattering matrix [Sr] due to 
a succession of duct sections as in (2.8) to replace the termination scat- 
tering matrix, [Tb], one essentially gets the scattering patterns due to a 
3-dimensional S-shaped duct just as for the 2-dimensional case described 
above. The solution as proposed in [2] can be used to find the scattered 


10 


field [E‘\ which is represented by a matrix equation as follows 

[£'] = {[«„] + [5.,][P]|Sr][P](|/]-[S !! ][P][5r][P])- , [S !1 ]}[B i ] 

( 2 . 10 ) 

with the scattering matrices defined as [S 2 i] coupling matrix from an 
external source to modal excitations, 

[Si 2 ] radiation matrix from modal fields to an external field point, 

[S 22 ] internal reflection matrix at the waveguide aperture, 

[P] the modal propagation phase delay matrix, 

[ E ’] incident field of the external source, and 

[5p] the effective reflection matrix at the end of the first section. 

All these matrices are well defined in [2] except the effective scattering 
matrix [Sr] which is defined in Equations (2.8) and (2.9). The transmission 
matrix, [Tjj] as appeared in (2.9), will be discussed in the next chapter. 
[E* ] and [f?‘] are column matrices with three elements representing the 
three vector components of E‘ and E l , respectively. Once [Sr) is found, 
Equation (2.10) can be used to simulate farfield patterns which will be the 
topic of Chapter 4. 
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Chapter 3 

Modal Transmission 
Coefficient for the Junction 
Between a Straight and an 
Annular Rectangular 
Waveguide Section 


To determine the junction transmission matrix, the modal fields of the 
straight and annular rectangular waveguides must be defined first; this is 
done in sections 3.1 and 3.2. The junction scattering matrix is then given 
in 3.3. 


3.1 Modal Field Structure of a Straight Rect- 
angular Waveguide 

Figure 3.1 shows the geometry of the straight rectangular duct section. The 
modal electric field E r± within this straight rectangular waveguide region 
( z < 0) may be represented in the usual manner [1,2, 3, 4] as: 

e t± = EE 

n m 

Here, e£ m and e^ m denote the transverse (to 2 ) electric vector mode func- 
tions for the TE nm and TM nm modes, respectively (where TE and TM 


A r A. e_ m e 

nm nm 




-f B r± (e 

1 nm V n 


1 * r 

~r ZP 
'nm znm 


) e 


TjPnr 


. (3.1) 
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Figure 3.1: Geometry of the straight rectangular waveguide section. 

are with respect to the i-direction). Likewise, the /3^ m and (3 r nrn denote 
the propagating constants of the TE nm and TM nm modes, respectively. It 
is noted that e T znm is the z-component of the TM nm modal electric field. 
The superscripts + and — in Equation (3.1) refer to modes propagating 
in the +z and — z directions, and the superscript r implies that the fields 
pertain to a straight rectangular waveguide section. It is also convenient to 
define the magnetic field H r± in the straight rectangular waveguide region 
following the representation for the electric field E T± in (3.1); thus 

H r± = EE 

n m 

(3.2) 

T| ie A’X, and B'„* in (3.1) and (3.2) are the modal coefficients which can 
be determined once the excitation is known. The ti^ m and zh T z ' nm are the 
magnetic vector mode functions for the TE nm modes, whereas h^ m is the 
magnetic vector mode function for the TM n m modes. Table 3.1 gives the 
explicit forms of the propagating mode functions in a straight rectangular 
waveguide. For the derivation and explicit expressions for these modes, one 


Kt (±C + iC.) e’*’--' + B (±K m ) 
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can refer to a well known text such as [5]. 


3.2 Modal Field Structure of an Annular Duct 

The modal fields for an annular duct also of rectangular cross-section, which 
is shown in Figure 3.2, can be established in a similar manner. Now, since 
the propagation is along the ±<£ directions, and 4> is not a constant vector, 
the fields E a * and H a± are initially defined as having components trans- 
verse to z and along z, respectively. Afterwards, these modal fields can be 
rearranged to give a more convenient set of modes; i.e., transverse to (j> and 
axial to <£, since ^ is the propagating direction. The modal electric field 
E a± and the magnetic field H a± can then be written as: 

= e e [<%* (« ± **,) s*** + K (n, ± **,) 

P 9 

(3.3) 

= EE [<# (±K + fat,) ***** + d£ {±K, + 

P 9 

(3.4) 

Again, e pq and e pq denote the transverse (to </>) electric vector mode 
functions for the TE pq and T modes ,and the h^ q and h° q denote the 
transverse (to <f > ) magnetic vector mode functions for the TE M and T M pqi 
respectively. The superscript a implies that the fields pertain in the annular 
duct section. The modal field structure of an annular duct will be discussed 
in more detail in the next section. 
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3.2.1 Modal Field Structure Calculations 

In an annular waveguide (as shown in Figure 3.2), there are always a finite 
number of propagating modes. These modes can be classified into sets of 
TM to z or TE to z, where the former set of modes has no H z field and 
the latter set of modes has no E z field. A combination of these two sets of 
modes gives a complete field structure inside the waveguide. In addition to 
these propagating modes, there are also an infinite number of evanescent 
modes. These evanescent modes decay exponentially along the waveguide 
axis. For the present analysis, the effects of these evanescent modes are 
ignored. 

TM (transverse to z) modes: 

The TM Z fields can be found via a z-directed magnetic vector potential, 
A = ^z. The magnetic vector potential satisfies the wave equation 

(V 2 + fc 2 )I = 0 (3.5) 
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where k is the free space wave number, 2tt/X. But, since z is a constant 
vector, the above equation can be reduced to a scalar form as: 


(V 2 + fc 2 )tf = 0. (3.6) 

This equation in cylindrical coordinates can be rewritten as: 

‘i d_ / d \ j_ f ' 

p dp \ dp) p 2 d<f> 2 dz 2 

The electric and magnetic fields are related to as [6]: 

’ p It (Pjt IT 10£ 

jk dpdz P p 6<f> 

' «* = -% (3-8) 

, -E, = |f (£ + *’) 4 H. = 0. 

The boundary conditions used to find 'f r are that all tangential electric field 
components are zero on the walls of the guide. In order to construct \t, let 
the separation of variables suggest a solution of the form: 


$ = 0. (3.7) 


<t(pAT) = R(p)iWZ(z). (3.9) 


Substituting Equation (3.9) into Equation (3.7) and introducing two sep- 
aration constants p x and u q , one will get three differential equations as 
follows: 


i a 2 

zdz iZ Pl 

(3.10) 

i d 2 $ 

— —v 

$ d<f> 2 9 

(3.11) 

-v]]R = 0. 

(3.12) 


The solutions to Equation (3.11) are 

$(<£) = e - ju '+ and 


(3.13) 


where modes propagate in the ±<f> directions. 
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The solution to Equation (3.10) is 

Z(z) = Acosp-iz + Bsinpxz. 
Since E<f, = the boundary condition 

E^{z = 0; z = /) = 0 


(3.14) 


(3.15) 


implies 


— (2 = 0; z = l) = 0. (3.16) 

Therefore, from (3.14) and (3.16), B = 0 and pi = e p. Hence, Z(z) can be 
rewritten as 


Z(z) — i4cos ^j-z, p = 0,1,*** (3.17) 

V 

Likewise, (3.12) is a Bessel’s equation and has solutions R(p ) which can be 
linear combinations of H^\k c p) and H^(k c p) with 

K = ^-(y) 2 (3*18) 

where H^\k c p) and H^(k c p) are Hankel functions of the first and second 
kinds, respectively, of order u q and argument k c p. From the boundary 
condition on the i directed E field, one obtains 


E z (p = a;p = b) = 0. 


(3.19) 




Next, one can utilize (3.19) and (3.8) to obtain 

dz 2 

zA 

i k 

Now R(p) can be constructed via (3.20) as 


k 2 + -=~z\ ¥ 


- (t)'] 


^(p = a; p = b) — 0. 


(3.20) 


R(p) = Yj [H^Kb)Hil\k,p) - Hl)Xk'b)HV(k'p)] . (3.21) 
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Note that R(b) = 0 has been employed. To satisfy ’® r (/> = a) = 0, the 
eigenvalues u g can be found numerically by solving the equation 


R(a) = 0. (3.22) 

Since R(p) in (3.21) depends on p and q, R(p) should be rewritten as R pq (p). 
Consequently, the eigenvalues u q found from (3.22) should be dependent on 
p and hence are renamed as Upg. 

Using linear combinations of these solutions and enforcing the boundary 
conditions gives the complete solution for the magnetic potential ^ which 
can be written as 

* = EEC^Wco S (^z)e--’‘--+ + C-JJ„(p)co S (^z)^ 

(3.23) 

where is related to the modal coefficients. 

PH 

The electric and magnetic fields inside the waveguide thus can be found 
from Equation (3.8). In addition, these modal fields can be put in a form 
to fit Equations (3.3) and (3.4). For example, the TM W electric modal field 
can be written as: 

2? = (<* ± (3.24) 


The subscript e means the field is TM or e type. Also by definition e pq are 
field components transverse to 0, thus from (3.8) one gets 


e* iv n* = zE e± 

•9 *P4 

sfL.e Tjv "* = 4>E£ 

4>pq ^ <Ppq 


(3.25) 


E l£r E m and E t>9 are found from (3.8) and (3.23) which are summarized 


in (3.28-3.30). 

The TM pq magnetic modal field Hf* can be put in a similar form, from 
Equation (3.4): 

H f* = (±h a pg + 4>hl pq )e^^. (3.26) 
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where 



r *;,,«** = pn% 

l = 4 .H&. 

( 3 . 27 ) 

The individual electric and magnetic field components can be 

summa- 

rized as follows: 



J 7 »e± 

^ppq 

Apq T (?) sin (? 2 ) R 'M e * iUr '* 

( 3 . 28 ) 

jTiei 

^H>q - 


( 3 . 29 ) 

r»e± 

&zpq — 

~ a «T [*’ - (f)] cos (t z ) 

( 3 . 30 ) 

I TC± 

11 

ppq 


( 3 . 31 ) 

rre± 

U 4>pq ~ 

Apqj cos R! pq {p)e* 3 ''*' 4 ’ 

( 3 . 32 ) 

Hit, = 

0 . 

( 3 . 33 ) 


The normalization constant Apq is found by normalizing the power car- 
ried by the pq th mode to unity. This is done by integrating the Poynting 
vector over a cross section of the guide and setting this equal to unity, 


/ X te x *») ■ * ds = L 

Substituting (3.25) and (3.27) into (3.34), A pq can be found as 


where 


. = {? 

1 1 ,p 

TE (transverse to z) modes: 


= 0 
/ 0 . 


(3.34) 


^ (x35) 


(3.36) 


The derivation for the TEpq fields is very similar to that for the TM pq fields. 
They can be found as follows: 
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Let F = be the electric vector potential satisfying the wave equation 


(V 2 + fc 2 ) F = 0 


(3.37) 


which can be reduced to a scalar form as in Equation (3.6). The electric 
and magnetic fields are related to this potential as [6] 

( j? it _ i dj* 

&P — pd<t> Ij p — jkZo dpBz 

* — jizl \ d&k (3.38) 

. e, = o a, = (& + *>) *. 

Again, using the boundary conditions that electric field components 
are zero on the walls, these field components can be put in a form as in 
Equations (3.3) and (3.4), such that 

Ef = « ± hi,) e****, (3.39) 


and from (3.38) one gets: 



= 

= 


(3.40) 


The subscript h in (3.39) and the same superscript in (3.40) means the field 
is TE or h type. Furthermore, E fj* and E ^ are summarized in (3.43) and 
(3.44), respectively. It is similar for the TE M case 



(±h a p ' g + 4>h a + pq ) 




(3.41) 


where 


{ 


= zH + pff “ 




'PPQ 






zpq 

h± 

4>pq' 


ppq 


(3.42) 


The individual field components can be summarized as follows: 

£& = B„^sin(^ 2 )fi M ( ( ,)«W''-* (3.43) 
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where 


j7»h± 

- 

-Bpqjs\n{~z 


(3-44) 

77l/l± 

“ 

0 



(3.45) 

rrh± 

**PP<1 “ 

-•-AC?) 

-(?•: 


(3.46) 

rr^± _ 
U 4>P<1 ~ 

*-A,(7 )■ 

108 (7 ! ) 

R.,(c)^ , "'' , 

(3.47) 

rr/i± 

**zpq ” 

; -VU* 2 - 

(^)] sin (*r*) 

(3.48) 

Rpq{p) 



(3.49) 


■ iRW 


K = 

The normalizing constant, J is 

l 


B pq — 


2fcZ 0 M 


with R P q{p) defined in (3.49). 


(3.50) 


-(t) 1 (3 - 5i) 
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3.3 Junction Transmission Matrix 


In this section, the junction transmission matrix from a rectangular waveg- 
uide to an annular waveguide is discussed. When the modal fields are 
transmitted from a rectangular duct section to an annular duct section, the 
field representation must be transformed from rectangular modes to annu- 
lar modes. Since both of these modal field structures are well-defined, the 
question left is to find the suitable excitation coefficients for the transmit- 
ted annular modes. A junction transmission matrix performs this function; 
it finds the annular transmitted fields once the incident rectangular modal 
fields are known. This transmission matrix is obtained by employing the 
incident modal fields to define the equivalent surface current at the junc- 
tion aperture. The latter equivalent sources then generate the transmitted 


annular fields. 

It is noted that the modal coefficients A and for the straight 
rectangular guide as in (3.1-2), and together with for the annular 
guide as in (3.3-4) can be represented as column vectors by 


TE Rectangular Guide Modes: 

TM Rectangular Guide Modes: 

TE Annular Guide Modes: [C^] 

TM Annular Guide Modes: [D£±] 


a t ± 


Ar± 

a nm j 

Bit 


B-Jm 
ctt ' 


C a± 

PQ 

D a± 


D a± 

U PQ 


(3.52) 

(3.53) 

(3.54) 


(3.55) 


If the rectangular modal fields are incident on a rectangular /annular duct 
junction, the modal coefficients for the transmitted fields in the annular re- 
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gion are related to the incident modal coefficients by a transmission matrix, 
T\ 2 , as 


' tel ' 


'[rJU] te„ m ] ' 


' ' 

. te + ] . 




. 1%] . 


Pial 


(3.56) 


where the subscript of [T 12 ] means propagation is from duct section 1 to 
duct section 2. In this case, section 1 is a straight rectangular section and 
section 2 is an annular duct section. Futhermore, the meaning of r^. nm is 
the following. A rectangular TM nm (or e type) mode with modal amplitude 
B which is incident at the junction is partly transformed into an annular 
TE„ (or h type) mode with modal amplitude as 


= 


phe D r + 
A pgjnm^nm 


(3.57) 


and the rest is transformed into an annular TMpq mode with modal co- 
efficient The + sign in the above equation defines the propagating 

direction towards the termination as indicated in Figure 1.1. When prop- 
agation reverses direction, the modal fields from an annular duct would be 
incident on the junction. The new transmission matrix [T 21 ] can be shown 
by reciprocity to be the transpose of the transmission matrix [T i2 ]. This 
relationship can be written as: 


Pii] = [r„]* 


(3.58) 


where t is a transpose operator. 

The junction transmission matrix [T , I2 ] as shown in Equation (3.56) has 
four submatrices. Each of these submatrices carries a different function. 
For example, the submatrix [r^. nm ] transforms the incident TM nm rectan- 
gular modes into the transmitted TE pq annular modes. These submatrix 
elements can be found by employing the incident modal fields to define 
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equivalent incident electric and magnetic surface currents in the aperture 
formed by the junction via the Kirchhoff approximation and then evaluat- 
ing the amplitudes of the transmitted modes which are produced by those 
equivalent currents as follows [1,2]: 

k = 4 [/ L • f - s * ■ ") H (3 - 59) 

where B°I is the transmitted modal excitation which is (7“ + or Dtf for 

PH PH Ph 

TE pq or TM pq transmitted, respectively. and H £+ are the transmitted 

modal fields which can be represented by Equations (3.28) to (3.33) for 
TMpq transmitted and by Equations (3.43) to (3.48) for TEpq transmitted. 
These fields can be expressed as 
for TM pq transmitted: 





4*7 + 


iElt, + + ■t’K, 

(3.60) 

K, + 



(3.61) 


e°' -(- ie'jL. 



(3.62) 

K + M* 


**% + ?*& + **& 

(3.63) 


Note that the phase factor has been suppressed because its effect 

will be included later in the propagating modal phase path delay matrix 
[P]. Furthermore the expressions for the annular modal fields components 
Etpq, E^pq, Efpq, H^, H ^ and H£ pq are summarized in Equations (3.28) 
to (3.33) for the TM case and in Equations (3.43) to (3.48) for the TE 
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I 

I 


case. J and M in (3.59) are the equivalent incident electric and magnetic 
currents, respectively. They can be written as 

J = hxH (3.64) 

M = Exn (3.65) 


where h is the unit aperture normal vector pointing in the direction of 
propagation of the incident modal field. The E and H in (3.64) and (3.65) 

—t — * 

are found via the Kirchhoff method as mentioned earlier; hence, E and H 
are assumed to be given by the TM and TE incident fields as follows: 

(1) for TM nm incident: 


E 

H 


^ nm znm 


~Km 

nm 


and (2) for TE nm incident: 


E 

H 


c nm 


~h r ' m + zKi 


(3.66) 

(3.67) 


(3.68) 

(3.69) 


The field components for both the TE nm and TM nm cases are summarized 
in Table 3.1. 

Each of the four submatrices can be represented by a combination of 
two sets of Equations. The first set of equations is chosen either from (3.60) 
and (3.61) for the transmitted annular TM pq fields or from (3.62) and (3.63) 
for the transmitted annular TEpq fields. The second set of equations com- 
prises either of (3.66) and (3.67) for the incident rectangular TM nm mode 
or of (3.68) and (3.69) for the incident rectangular TE nm mode. For ex- 
ample, substituting the first set of Equations (3.60), (3.61) and the second 
set (3.68), (3.69) into Equation (3.59), one gets the transmission submatrix 
r;$. nm which transforms a rectangular TE nm mode into an annular TM m 
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mode. The four possible combinations of these two sets of equations gener- 
ate the four submatrices in the transmission matrix. The formulations for 
these submatrices are summarized in Table 3.2. 

There is one junction transmission matrix for each junction in the duct 
structure as shown in Figure 1.1. However some of the junction matrices 
are equal due to symmetry. These relationships are as follows: 

[T„] = [T 32 ] = [T 23 ] f (3.70) 

and, if r 2 = r 4 in Figure 1.1, then 

Pnl = [r„] = [T„]*. (3.71) 

There is a reversal of coordinate systems at the junctions between sec- 
tions 3 and 4. Therefore, [T 4B ] ^ [T 43 ]. However, it can be shown that there 
is a simple relationship between the elements of the matrices [Xi 2 ] and [jT 34 ]. 
For both the rectangular TE and TM cases, this relationship is given by 

[T 34 ] nm = — (— l) n [T 12 ] nm (3.72) 

where nm is the rectangular mode number in the width (^-direction) and 
in the height (y-direction) dimensions, respectively. Due to the symmetry 
relations of the transmission matrices, if the inner radius of the curved 
sections are equal, i.e., r 2 = r 4 , then only [Ti 2 ] is needed to be evaluated in 
order to find [j 4] in (2.9). 


I 
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Table 3.1: Properties of Modes in Rectangular Waveguide 


K 

h x 

ky 


TE„ 


y nz. -f m . 

nm ia cos n a x cos m b y 

JPnm 


= 0 


= N nm Y nm n a sin n a x cos m b y 
= N nm Y nm m. b cos n a x sin m b y 
= n nm m b cos n a x sin m b y 
— —N nm n a sin n a x cos m b y 


TM nm 

0 

nr ni+m? 

-N nm sin n a x sin m b y 

-N nm Y nm m b sin n a x cos m b y 
N nm Y nm n a cos n a x sin m b y 
N nm Tl a COS Tl a X S1TI TTl b y 
N nm m b sin n a x cos m b y 


Y nm = Y 0 /3 ntn /k 


Y 0 k/f3 nm 

o < x < a, 0 < y < b 


Note that for both modes: 
W nm 4 [2e <m c <mt yi xm a6(n a + m^)] 2 

/?nm = [k 2 -n\- m \\ 2 

n a = Tin /a 


ton 


m b = mn/b 


2 n = 0 
1 n ^ 0 
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Table 3.2: Transmission Matrix [T 12 ] from a Rectangular Duct to an An- 
nular Duct 

(1) TM nrn (rectangular) modes in, TM M (annular) modes out 

| %A„ (k 2 - (f ) 2 j N nm Y nm m b / r r ; sin n a xR pq {p)dp 

|f A M (&■) N nm Y n mn a cos n a xR' pq {p)dp 
{ - \v Pq N n mTn b A pg l sin n a xR pq { P )d P 

if p = m 

and =0 if p ^ m 
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Table 3.2 continues 


(2) TM nm (rectangular) modes in, TEpq (annular) modes out 



-;5fc ( fc2 ” (“)*) N ™ n ‘fr\ a cos n a xR pq {p)dp 

-{§? (*r) ^ nmrrib /r» B sin n a xR' pq (p)dp 
- k 4 i^N nm Y nm n a B pq / r r ; J cos n a xR pq (p)dp 


if p = m 

and = 0 


if p m 


with 


7T — 1,2,3,.** 
m = 1,2,3, .. . 


and 


p — 1,2,3,... 
9 = 1, 2, 3,... 


B pq are defined in Equation (3.51). R^R'^ and v m are defined in the 
section oiTE annular modes. 
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Table 3.2 continues 


(3) TE nm (rectangular) modes in, TM pq (annular) modes out 



| (k 2 - N nm Y nm n a / r r * sinn a xR pq (p)dp 

< \eApg N nm Y nm m b cos n a xR' pq (p)dp 
\eA pq u pq N nm n a / r r “ J sin n a xR pq {p)dp 


if p = m 

and =0 if p ^ m 


with < 


n = 0, 1, 2, 3, . . . 
m = 0,1, 2, 3,... 
n = m — 0 excepted 
p = 0, 1, 2, 3, . . . 

^ — 1,2,3,... 


and e 


f 2, P = 
1, P = 


= m = 0 

m ^ 0. 


ApgiRpg,^ and t/p, are defined in the section of TM annular modes. 
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Table 3.2 continues 


(4) TE nm (rectangular) modes in, TE n (annular) modes out 

“IS* ( fc2 ” (*r) 2 ) N n m m b l ;• cos n a xR pq (p)dp 
[ r M ;nm] = < - J g* (*= ) N nm n a / r p ; sin n a xR' M (p)dp 

k 4 ^pg^nm^'nm'^'b^pq frf, ~ COS Tl a xRpq(p)dp 

if p = m 

and =0 if p ^ m 

n = 0, 1,2, 3,. . . 
m = 0,1, 2,3,... 
with < n = m = 0 excepted 
p = 1,2,3,... 
k 9 = 1,2,3,... 

B pq are defined in Equation (3.51). i2 pg ,i2^ and t'p, are defined in the 
section of TE annular modes. 
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Chapter 4 

Results and Discussion 


In this section the fields backscattered from a rectangular, S-shaped duct 
will be presented in the form of backscatter vs. aspect angle plots. The 
geometries of an S-shaped duct and a straight duct with rectangular cross- 
sections are illustrated in Figures 1.1 and 3.1, respectively. The scattering 
effects included in this analysis are the internal reflection at the duct cavity 
termination and the external rim diffraction at the duct opening. The 
structural scattering due the exterior duct surface is not included. The 
emphasis of this section will be focused on comparing the backscattered 
fields of an S-shaped duct with those of a straight duct. It is expected 
that the backscattered patterns will be different because the curvature of 
an S-shaped duct redistributes the coupled energy in a way that differs 
from a straight duct. Since both the straight and S-shaped ducts in this 
analysis are perfectly conducting, none of the energy coupled into the duct 
is lost. Although the backscattered pattern will change depending on the 
duct geometry, the average dB level will remain about the same. 

First, as a check to the analysis of an S-shaped duct, the backscattered 
patterns of a 5 section S-shaped duct is obtained in the x — z plane. As 
shown in Figure 4.1, the lengths of the last 4 sections of this S-shaped duct 
approach zero. In other words, the geometry of this S-shaped duct is no 
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different than a straight duct. The height of this duct is 6 inches and the 
width is 4 inches. At 10 GHz, this translates to 99 propagating modes in the 
rectangular section and 105 modes in the annular section. The backscat- 
tered E-theta plot and the E-plii plot are shown in Figures 4.1 and 4.2, 
respectively. The corresponding patterns from a straight duct are shown in 
Figures 4.3 and 4.4. Comparing the backscattered .E-theta patterns of Fig- 
ures 4.1 and 4.3 and comparing the backscattered E-phi patterns of Figures 
4.2 and 4.4, one finds that they are virtually the same. Consequently, this 
analysis successfully reduces a more complicated S-shaped duct geometry 
into a straight duct geometry as far as backscatter patterns are concerned. 

Figures 4.5 and 4.6 show the backscattered patterns of an S-shaped duct 
in the x — z plane. Figure 4.5 also defines the geometry of this S-shaped 
duct. Comparing these two patterns with those due to a straight duct as 
shown in Figures 4.3 and 4.4, one finds that the S-shaped duct changes the 
pattern shape. The coupled energy is redistributed after it radiates back 
out although the average dB level remains about the same. Moreover, in the 
x — z plane, the S-shaped duct does not have symmetry in geometry as the 
straight duct does. Thus the backscattered patterns of the S-shaped duct 
in this plane are not symmetric, contrary to the patterns due to a straight 
duct; these effects can be seen from Figures 4.3 to 4.7. The backscattered E- 
theta field of this S-shaped duct, in Figure 4.5, also compares very well with 
the TM backscatter of a 2-dimensional S-shaped duct as shown in Figure 
4.7. This is expected because only the pattern in the principal (horizontal) 
plane of this 3-dimensional duct is considered. The field patterns due to 
the other principal (vertical or y — z) plane is shown in Figure 4.8. It can be 
seen that in this plane, the pattern is symmetric because of the symmetry 
of the structure. 

Figures 4.9 to 4.11 show the backscattered patterns of the S-shaped duct 
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due to different conical pattern cuts. The theta angle is fixed at 10,30 and 
45 degrees for Figures 4.9, 4.10 and 4.11 respectively. The phi angle changes 
from 0 to 360 degrees for each of these figures. It can be seen that the dB 
levels of these plots are about the same. However, the larger the angle 
theta, the more the pattern fluctuates. Such characteristics can also be 
seen in Figure 4.6 which shows nine conical pattern cuts at every 5 degrees 
from 5 to 45 degrees. Again, this is expected because a small theta angle 
leads to a small conical pattern variation in circumference. For example, 
in the extreme case, if theta is zero, the whole conical pattern region will 
be reduced to a single direction. As a result, a conical pattern with a small 
theta will not be influenced as much by the radiation from many dominant 
higher order modes as a conical pattern with a larger theta. 
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Figure 4.7: Geometry and TM backscatter plot of a 2-dimensional S-shaped 
duct at 10 Ghz. 
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Backscattercd E-Theta Plot 

Figure 4.11: Conical pattern cut at 6 = 45°. 
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